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SUMMARY 

A  theory  of  sound  radiation  from  an  infinite  cylindrical  shell  which  is  stiffened  periodically  by  ribs 
is  given.  Both  shell  and  ribs  are  laminated  composites.  The  shell  is  modelled  by  the  differential 
equations  of  a  thin  shell  theory.  The  flexible  ribs,  which  exert  both  forces  and  moments  on  the 
shell’s  surface,  are  modelled  by  thin  conical  shell  finite  elements;  degrees  of  freedom  other  than  those 
at  the  shell’s  surface  are  removed  by  receptance  methods.  Excitations  are  time-harmonic  mechanical 
forces  at  arbitrary  points  on  the  shell’s  surface  and  ribs,  and  acoustic  monopoles  in  the  exterior  and 
interior  fluids. 
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SOUND  RADIATION  FROM  INFINITE  LAMINATED  COMPOSITE  SHELL  WITH  PERIODIC 
RIB  STIFFENING 

1.  INTRODUCTION 

An  increasing  use  of  laminated  fibre  reinforced  materials  as  components  in  marine  vehicles,  see,  for 
example.  Smith  [1],  means  that  well-established  theoretical  models  based  on  isotropic  materials  must 
be  re-worked  to  a  certain  extent  to  include  anisotropy.  The  particular  problem  studied  here  is  the 
sound  radiation  from  a  laminated  composite  cylindrical  shell  which  is  stiffened  periodically  by 
identical  ribs  with  flexible  cross-sections.  The  aim  is  to  extend  previous  theoretical  formulations  to 
cover  laminated  composite  materials,  flexible  ribs  of  arbitrary  cross-section,  interior  and  exterior 
acoustic  excitation,  and  mechanical  excitation  on  a  rib. 

Sound  radiation  from  stiffened  isotropic  thin  shells  has  been  studied  theoretically  by  several  authors. 
Con^’ersion  of  subsonic  non-radiating  wavenumbers,  at  the  ribs,  to  supersonic  radiating  wavenumbers 
is  recognised  as  an  important  source  of  sound  radiation  with  quasi-resonance  characteristics.  In 
Bemblit  [2]  the  rib  is  modelled  as  a  simple  circular  beam  which  interacts  with  the  shell  through  radial 
forces  only;  excitations  are  radial  forces  on  the  shell’s  surface.  In  James  [3]  a  matrix  formulation  is 
given  which  allows  the  rib  (and  excitations)  to  exert  axial,  tangential  and  radial  forces  together  with 
a  meridional  moment  on  the  shell’s  surface;  the  rib  is  modelled  as  a  circular  beam  whose  centroid 
may  be  offset  from  the  shell’s  surface;  subsequent  unpublished  work  modelled  a  flexible  rib  by  the 
finite  element  method.  In  Burroughs  [4]  a  doubly  periodic  arrangement  of  ribs  is  considered;  the 
ribs,  which  are  assumed  to  interact  with  the  shell  through  radial  forces  only,  are  modelled  as  beams 
with  coupled  rad'al  and  tangential  motion;  excitations  are  radial  forces  on  the  shell’s  surface.  In 
Burroughs  and  Hallender  [5]  the  previous  work  has  been  extended  to  cover  several  different  types  of 
mechanical  excitation  on  the  shell’s  surface.  In  Hodges  et  al.  [6]  free  wave  propagation  on  an 
unloaded  cylinder  is  examined;  in  order  to  give  the  ribs  some  measure  of  flexibility,  the  web  of  their 
"T-rib”  is  modelled  with  cubic  variation  for  its  axial  displacement;  this  model  has  been  substantially 
updated  by  Brazier-Smith  and  Scott,  unpublished  work  1988,  Topexpress  Ltd,  to  include  far  field 
sound  radiation  and  certain  hydroacoustic  noise  sources  as  excitation.  References  to  other  relevant 
work,  .  uticularly  on  wave  propagation  in  and  sound  radiation  by  stiffened  plates,  are  contained  in 
the  articles  cited. 

The  geometry  is  shown  in  Figure  1.  A  thin  cylindrical  shell  is  stiffened  by  periodically  spaced 
identical  axisymmetric  ribs.  The  cross-section  of  a  rib,  shown  as  a  "T-rib”  in  the  diagram,  is 
composed  of  thin  elements.  Excitations  are  time-harmonic  mechanical  forces  on  the  shell’s  surface 
and  ribs,  and  time-harmonic  acoustic  monopoles  in  the  exterior^  and  interior  acoustic  fluids.  The 
effect  of  sound  scattering  by  the  finite  sized  ribs  is  neglected,  ie  for  the  purposes  of  calculating  the 
acoustic  pressure  inside  the  shell  the  ribs  are  regarded  as  line  attachments.  Both  shell  and  ribs  are 
laminated  composites  comprising  fibre  reinforced  thin  layers  with  arbitrary  sta'^king  angles.  In 
Figure  2(a)  is  shown  the  through  thickness  lamination,  as  a  planar  stack  of  layers;  in  Figure  2(b)  is 
shown  the  stacking  angle  of  a  typical  composite  layer,  the  directioni  {<l>,z)  on  the  shell’s  surface  being 
shown  in  parentheses. 

The  work  described  herein  was  done  under  item  AS02EH14  of  the' Strategic  Research  Programme. 
It  contributes  to  the  first  milestone.  1 
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2.  MATHEMATICS  OF  PROBLEM 

2.1.  Fourier  Transforms 

Field  quantities  are  represented  by  Fourier  integral  transforms  in  the  axial  direction  and  Fourier  series 
transforms  in  the  circumferential  direction;  for  example,  the  function  P(r,^,z)  is  expressed  as 

-  (l/2n)  E  /"  P{r^,a)e'"’da,  (2.1.1) 

where  the  spectral  quantity  is  given  by  the  inverse  transform 

Pir,n,a)  =  (l/2it)j‘%-‘^|‘ (21-2) 

Fourier  transforms  facilitate  solution  of  linear  differential  equations  by  replacing  them  with  linear 

algebraic  equations.  In  these  and  subsequent  equations  the  time  variation  e'**"  is  suppressed. 
Hereafter,  bars  on  spectral  quantities  will  be  omitted  to  simplify  notation;  the  arguments  of  a  function 
are  sufficient  to  identify  it  as  a  Fourier  series  or  Fourier  integral  transform. 

2.2.  Spectral  Response  of  Shell 

Let  u(^,z)«{ttj(0,2),a^(0,2),B,(0,2),^r^(^,?)y  denote  a  4  x  1  column  vector  of  axial,  tangential  and 
radial  displacements  and  meridional  rotation  of  the  shell’s  midsurface.  From  now  on,  references  to 
a  displacement  vector  implicitly  includes  the  meridional  rotation,  and  references  to  an  excitation 
vector  implicitly  includes  the  meridional  moment.  The  spectral  displacements,  using  general  thin  shell 
and  rib  equations  of  motion,  have  been  obtained  by  James  [3]  as  the  4  x  4  matrix  equation 

u(/i,a)  ■  D(rt,«)E(n,a)-(l/<i)D(«,o)B(rt)  x 

Uim'£h:n.«*27tql<mn) 

.1  4—* 

where  u(R,a)  ■  is  the  4  x  1  column  vector  of  spectral  displacements 

of  the  shell’s  midsurface;  D(n,a)  is  a  4  x  4  spectral  receptance  matrix  whose  elements,  given  in 
Annex  A,  depend  only  on  the  shell  and  fluid  parameters;  E(R,a)  is  a  4  x  1  spectral  excitation  matrix 
whose  elements  are  obtained  from  the  external  excitation  matrix 

in  which  EJiip^)  and  £,(^,z)  are  axial,  tangential 

and  radial  external  stresses  and  Af^(^,z)  is  the  meriodional  moment  excitation  per  unit  area;  d  is  the 

distance  between  ribs;  B(n)  is  a  4  x  4  spectra!  stiffness  matrix  of  a  rib  whose  elements  when  it  is 
modelled  as  a  beam  are  given  in  reference  [3]  and  whose  elements  when  it  is  modelled  by  thin 
axisymmetric  conical  shell  finite  elements  are  derived  according  to  principles  described  in  Annex  B; 
I  is  a  4  X  4  identity  matrix.  Inspection  of  the  above  equation  shows  that  the  spectral  displacement 
matrix  is  composed  of  two  terms:  the  first  term,  of  axial  wavenumber  o  is  the  .spectral  displacement 
of  the  shell  in  the  absence  of  ribs;  the  second  term  comprises  discrete  axial  wavenumbers  a  ■^2xq/J 
due  to  various  interactions  between  the  periodically  spaced  ribs,  via  the  shell  and  its  interior  and 
exterior  fluids. 


i-t 


5^  D(,n,a+2rtqld)E(n,a-*-2itqld) 


«— 


(2.2.1) 
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2.3.  Far  Field  Sound  Radiation 


The  sound  radiation  in  the  exterior  fluid,  due  to  a  prescribed  spectral  radial  displacement  u/n,a)  of 
the  shell’s  surface, 

and  its  stationary  phase  value  in  the  acoustic  far  field,  /{-«•, 


Pf{R,Q.<P)  =  (2-3.2) 


are  results  that  are  a  standard  to  acoustics  (see,  for  example,  references  [7]  and  [8]).  In  these 

equations  is  the  density  of  the  exterior  fluid;  is  the  acoustic  wavenumber 

in  the  exterior  fluid  where  c,  is  sound  speed;  /f,  is  the  Hankel  function  and  the  prime  on  the 
Hankel  function  denotes  differentiation  with  respect  to  its  argument;  (R,d,4>)  are  spherical  coordinates 
related  to  cylindrical  coordinates  by  3:=Rcos(^,  r=Rsm{d):  ujinjt^cosd)  is  the  third  element  of  the 
spectral  displacement  vector  of  equation  (2.2.1)  evaluated  at  wavenumber  a=Jt^cos5.  Provided  that 
there  is  some  dissipation,  however  small,  in  the  shell  and  ribs,  the  method  of  stationary  phase  is  valid 
formally  for  all  values  of  6  and  4>- 

2.4.  Mechanical  Excitation  on  Shell  Surface 


The  basic  external  excitation  t tress  on  the  shell’s  surface  is  the  point  force  and  meridional  point 
moment  excitation  vector,  F>«{F,/'^,,Af^,  located  at  coordinates  (fl,^,,z^.  In  this  case 

E(0,z)  -  ¥i(z-z^)ii<p  -<f>^)la,  (2.4.1) 


and  the  spectral  excitation  is  obtained  from  equation  (2.1.2)  as 

E(n,a)  -  Pexp(-in^^-iaz,)/2i:a.  (2.4.2) 


2.S.  Mechanical  Exciution  on  Rib 

A  rib  is  modelled  by  conical  shell  finite  elements,  as  shown  in  Figure  3,  the  theory  being  given  in 
Annex  B.  If  there  are  nodes  on  a  rib  cross-seaion  there  are  4  Af  degrees  of  freedom,  which  could 
considerably  increase  the  size  of  the  matrices  of  equation  (2.2. 1)  unless  unwanted  degrees  of  freedom 
are  first  eliminated.  It  is  assumed  that  the  rib  is  connected  to  the  shell,  along  a  circumference,  at 
node  p,  say,  where  there  are  four  displacements  and  four  excitations  which  match  the  shell 
displacements  and  exciutions.  If  a  mechanical  point  force  F  is  applied  at  node  9  on  a  flexible  rib, 
then  Annex  B  shows  how  the  method  of  receptance  coupling  enables  it  to  be  replaced  by  an  equivalent 
shell  excitation  at  the  point  p  where  the  rib  connects  to  the  shell.  This  shell  excitation  is 
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E(n,a)  -  R^(n)*‘Rj^(n)Fexp(-in^,-ioz,)/2jia,  (2.5.1) 

where  R^(n)  and  R^(n)  are  4  x  4  rib  receptance  matrices.  The  columns  of  R^(n)  are  the  four 
Fourier  series  harmonic  displacements  of  the  rib  at  p  due  to  unit  harmonic  excitation  at  each  in  turn 
of  the  four  excitation  degrees  of  freedom  at  q .  When  the  points  p  and  q  are  coincident  it  is  evident 
that  equation  (2.S.1)  reduces  to  equation  (2.4.2).  Thus,  provided  the  excitation  is  applied  at  a  single 
node  on  a  flexible  rib,  the  method  of  receptance  coupling  allows  reduction  of  the  number  of  rib 
coupling  degrees  of  freedom  from  4  V  to  four. 

2.6.  Acoustical  Excitation 

A  point  source  located  in  the  interior  fluid  at  coordinates  has  free  field  pressure 

PjirAz)  •  pjtxf(ik^^)IR^  (2.6.1) 

where  -  (z-z,)*+r*+r*-2rr^cos(^-^^;  >b>/C|  is  the  wavenumber  and  e,  is  sound  speed. 

Elements  of  the  excitation  vector  E(n,a)  are  zero  except  for  [7] 

■  2p,[j|,|(Y^,)/Y|a^|,|(Y,a)jcxp(-in^,-f«z,),  (2.6.2) 

where  i 
a.3.2). 

When  the  source  is  located  in  the  exterior  fluid  at  coordinates  the  elements  of  the  excitation 

vector  E(n,a)  are  again  zero  except  for  [7). 

£,  -  2p,[ff,,,(Y/,)/Y,a^r|',|(Y,a)]exp(-w^,-ioz.).  (2-6.3) 

The  far  field  sound  radiation  due  to  the  shell’s  response  is  also  given  by  equation  (2.3.2),  which  must 
now  be  augmented  by  the  far  field  pressure  of  the  source 

•  (p,/R)exp[tf:^(R-z,cos0-r^sinft:os(^-^,))j,  (2.6.4) 

and  the  far  field  pressure  of  its  rigid  boundary  reflection  [7] 

P,f(K0,4>)  -  -(p,//0exp[a:,(R-z,co6^]  x 

(-0'"'[*^pi|(*/»in^/f|,,(*:/^“^//f|',|(t,flsin^]exp(irt(^-^,)).  (2.6.5) 


*-«*).  The  far  field  sound  radiation  due  to  the  shell  response  is  given  by  equation 
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NUMERICAL  EXAMPLES 


■j 

3.1.  General 

A  Fortran  program  has  been  written  to  calculate  the  far  field  sound  radiation  from  a  ribbed  laminated 
shell  for  both  acoustical  and  mechanical  time-harmonic  excitation.  In  order  to  provide  a  test  of  the 
capability  of  this  program  an  air,  pj-1.21,  Cj«343,  filled  cylinder,  radius  a  -  0.10,  thickness 
h  -  0.002 ,  excited  by  a  unit  rms  mechanical  point  force  and  radiating  into  water,  p^-1000..  c,-1500, 
has  been  considered.  The  cylinder  material  is  either  isotropic  steel, 
E  •  195,  G  •  75.6,  v-0.29,  fi*nO0,  or  carbon  fibre  reinforced  plastic  (CFRP), 
£^-202,  £^“7.4,  G^"2.74,  v^-0.320,  p-1520  when  reinforced  in  the  circumferential  direction 

and  an  interchange  of  ^  and  z  is  required  when  reinforced  in  the  axial  direction.  The  ’T-rib” 
dimensions  are  web  thickness  0.002,  web  length  0.09,  flange  length  0.08  and  flange  thickness  0.002; 
the  rib  spacing  is  <f  >  0.04  and  the  rib  material  is  either  steel  or  circumferentially  reinforced  CFRP. 
Ail  material  constants  are  in  SI  units  with  Young’s  and  shear  moduli  in  GN/m‘.  Dissipation  in  the 
cylinder  and  ribs  is  allowed  for  by  setting  Young’s  and  shear  moduli  to  complex  values,  viz 
£«£(l-iH,)  and  G«G(l-itip,  etc  where  the  hysteretic  loss  factor  and  are  all  chosen  as  0.02 
for  the  numerical  examples. 

The  angular  variation  of  the  far  field  sound  radiation  exhibits  a  complicated  variation  with  the  angles  d 
and  ^  due  to  the  dispersive  nature  of  wave  propagation  in  ribbed  cylinders.  A  feature  of  this 
dispersion  is  that  frequency  locations  of  peaks  in  the  radiated  sound  spectra  move  a.s  the  angle  of 
observation  changes.  In  order  'o  smooth  out  this  complicated  spatial  variation  the  acoustic  power  is 
plotted  in  Figures  4-5.  The  far  field  pressure  is  always  of  the  form 

• 

Pj  {R,e,<P)  -  exp  ]C  “P  (3.1-1) 

•••• 

and  it  is  not  difficult  to  show  that  the  acoustic  power  is 

i*,- (2n/p.c.)/’{E 

where  a  customary  multiplying  factor  of  1/2  has  been  omitted  as  it  is  convenient  to  think  of  the 
excitation  as  rms.  This  equation  is  integrated  numerically  by  an  extended  Simpson’s  rule.  The  power 
level,  in  dB  reference  1  picowatt,  is  defined  as  10Iog,o(Pp+ 120.0.  For  computations  of  acoustic 
power  a  maximum  of  10  circumferential  harmonics  (n),  12  axial  harmonics  (q)  and  361  theta 
increments  were  used.  Symmetry  in  the  ranges  6=(0,90)  and  fl  =  (90,180)  degrees  was  used  to 
reduce  computational  times. 

3.2.  Isotropic  Cylinder 

In  Figure  4(a)  is  shown  sound  power  levels  of  an  unribbed  steel  cylinder.  When  the  excitation  is  a 
radial  force  ^e  spectrum  is  generally  smooth  and  it  can  be  demonstrated  that  it  is  substantially 
independent  of  the  shell  parameters.  When  the  excitation  is  an  axial  force,  power  levels  are  initially 
much  greater  than  those  of  tlie  radial  force,  and  thereafter  the  levels  are  within  4  dB  of  each  other. 
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A  detailed  physical  interpretation  of  features  in  the  plots  is  not  given  here,  but  it  is  worthwhile 
recording  that  certain  aspects  can  be  explains!  qualitatively  from  wavenumber  versus  frequency  plots, 
as  in  James  [9].  For  example,  such  plots  would  show  that  while  there  are  no  dominant  acoustically 
fast  radial  modes,  there  are  acoustic^ly  fast  axial  modes  cutting  on  at  0,  S,  10  and  IS  kHz  for  the 
It  <■  0,  1,  2  and  3  harmonics,  respectively.  Thus,  most  of  the  power  of  the  radial  force  coaies 
from  near  field  cylinder  distortion,  while  most  of  the  power  of  the  axial  force  from  comes  from 
supersonic  quasi-axial  waves  which  radiate  acoustic  energy  in  very  sharp  beams  centred  on 
coincidence  angles  which  vary  with  harmonic  number,  n. 

In  Figures  4(b)  and  4(c)  are  shown  sound  power  levels  of  the  ribbed  shell  when  the  excitations  are 
on  the  shell’s  surface,  at  a  rib  attachment  point,  and  on  one  of  the  ribs,  at  its  flange  centre  (node  4 
in  Figure  3),  respectively.  When  the  excitation  is  a  radial  force,  there  is  no  noticeable  difference 
between  the  plots,  as  the  frequency  is  too  low  for  web  longitudinal  resonances  or  flange  bending 
resonances  to  occur.  The  broad  hump  at  10  kHz  is  more  than  10  dB  above  the  corresponding  level 
of  the  unribbed  shell.  It  is  due  mainly  to  a  near  field  effect  caused  by  increased  circumferential 
stiffness  of  the  ribbed  shell,  together  with  conversion  of  subsonic  wavenumbers  at  the  ribs  into 
supersonic  rad  iating  wavenumbers,  the  latter  effect  being  more  noticeable  in  sound  pressure  versus 
frequency  plots  for  selected  d  and  ^  angles.  When  the  excitation  is  an  axial  force  there  is  a  large 
increase  in  power  level  when  the  excitation  is  applied  to  a  rib  web  rather  than  to  the  shell’s  surface, 
most  certainly  due  to  resonant  bending  of  the  web  crc’S-section. 

3.3.  Fibre  Reinforced  Cylinder 

In  Figure  5(a)  is  shown  sound  power  levels  when  the  shell  is  composed  of  axially  reinforced  CFRP. 
When  the  excitation  is  a  radial  force  it  can  be  demonstrated  that  the  power  levels  are  generally  only 
a  few  dB  above  those  that  would  occur  in  the  absence  of  the  shell.  When  the  excitation  is  an  axial 
force  power  levels  are  on  average  about  10  dB  less,  except  at  frequencies  below  5  kHz. 

In  Figures  5(b)  and  5(c)  are  shown  sound  power  levels  of  the  ribbed  shell,  the  shell  being  composed 
of  axially  reinforced  CFRP  and  the  ribs  of  circumferentially  reinforced  CFRP;  excitations  are  on  the 
shell’s  surface,  at  a  rib  attachment  point,  and  on  one  of  the  ribs,  at  its  flange  centre  (node  4  in 
Figure  3),  respectively.  When  the  excitation  is  a  radial  force,  there  is  no  noticeable  difference 
between  the  plots,  except  in  the  range  15*20  kHz  where  levels  are  higher  when  the  excitation  is  on 
the  rib;  due  presumably  to  the  onset  of  a  web  longitudinal  resonance.  The  broad  hump  at  10  kHz  is 
less  dominant  than  in  Figure  4.  When  the  excitation  is  ai.  axial  force  there  is  a  large  increase  in 
power  when  the  excitation  is  applied  to  a  rib  web  rather  than  to  the  shell’s  surface,  most  certainly  due 
to  resonant  bending  of  the  web  cross-section. 

Comparing  power  levels  in  Figures  4  and  5  shows  that  when  the  excitation  is  a  radial  force,  levels 
of  the  CFRP  shell  are  generally  just  above  those  of  the  steel  shell;  for  axial  force  excitation,  levels 
of  the  CFRP  are  generally  considerably  less. 

3.4.  Comparison  with  an  Exact  Theory 

It  is  generally  agreed  that  higher  order  shear  deformation  theories  are  necessary  for  accurate  static 
and  dynamic  response  calculations  of  laminated  composites,  which  are  usually  very  compliant  in 
transverse  shear.  Classical  plate  and  shell  theory  (as  used  herein)  are  thought  to  be  inadequate  for 
all  but  the  thinnest  of  composites;  thus,  the  numerical  results  shown  need  to  be  validated  by 
comparing  them  with  those  of  a  higher  order  theory.  Fortunately,  a  Fortran  program  [7]  is  availabie 
for  predicting  acoustic  radiation  from  cylindrically  laminated  composites  for  the  case  in  which  the 
layers  are  modelled  by  exact  equations  of  linear  anisotropic  elasticity. 
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For  the  isotropic  steel  cylinder,  exact  theory  give^  numerical  values  of  power  which  are  almost 
indistinguishable  from  those  shown  in  Figure  4(a).  For  the  axially  reinforced  CFRP  cylinder, 
differences  between  power  levels  of  exact  theory  and  those  shown  in  Figure  5(a)  are  hardly 
noticeable.  The  plots  of  power  obtained  from  exact  theory  are  not  reproduced  here.  As  the  thickness 
of  the  cylinder  is  increased,  without  other  constants  changing,  the  differences  between  exact  power 
levels  and  those  of  shell  theory  will  become  increasingly  apparent,  but  an  asses..ment  of  the  range  of 
applicability  of  classical  shell  theory  applied  to  laminated  composites  is  outside  the  scope  of  this 
theoretical  paper.  However,  provided  the  composite  is  sufficiently  thin,  it  is  surmised  that  classical 
shell  theory,  used  with  caution,  is  adequate  for  phenomenological  numerical  studies. 

The  receptance  matrix  S(n,a)  of  the  shell  and  the  stiffness  matrix  B(n)  of  the  rib  can  also  be  obtained 
from  the  exact  theory  described  in  references  [7]  and  [13],  respectively.  The  analysis  procedure 
would  not  be  difrlcult,  but  computer  times  would  be  an  order  of  magnitude  greater. 

4.  CONCLUDING  REMARKS 

Previous  work  on  sound  radiation  from  an  infinite  shell  with  periodic  stiffening  has  been  extended 
to  include  fibre  reinforced  materials,  flexible  ribs  of  arbitrary  thin-walled  cross-section,  mechanical 
excitation  on  a  rib  and  acoustic  excitation  by  interior  and  exterior  monopoles.  The  flexible  rib  which 
is  modelled  by  conical  shell  finite  elements,  usually  has  a  large  number  of  degrees  of  freedom;  the 
method  .  f  receptance  coupling  has  shown  how  unwanted  degrees  of  freedom  can  be  eliminated. 
Numerical  results  have  demonstrated  that  the  presence  of  flexible  ribs  or  anisotropy  can  substantially 
change  the  acoustic  power  spectrum  characteristics.  Lengthy  parametric  studies  which  include 
different  fibre  and  matrix  materials,  stacking  angles  and  stacking  sequences  will  be  required  to 
establish  physical  principles  of  practical  importance.  The  computer  programs  developed  as  part  of 
this  study  are  usefiil  tools  for  such  parametric  studies. 

For  numerical  investigations  in  which  the  shell  contains  a  finite  number  of  ribs  of  arbitrary  spacing 
and  differing  cross-sections,  the  method  of  dynamic  stiffness  coupling  as  used  by  James  flO]  to  solve 
a  scattering  problem  is  relevant.  For  analytical  investigations  of  the  physics  of  sound  radiation  and 
scattering  the  methods  of  Skelton  ( 1 1]  are  appropriate  for  simple  constraints.  For  finite  axisymmetric 
shells,  with  internal  axisymmetric  structure,  models  based  on  analytical  methods  such  as  in  Laulagnet 
and  Guyader  [12]  or  coupled  finite  element  and  Helmholtz  integral  equation  formulations  as  in 
James  [13]  enable  computations  of  low  frequency  sound  radiation. 
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FIGURE  1.  (a)  PERIODICALLY  RIBBED  SHELL;  (b)  COORDINATE  SYSTEMS; 
(c)  POSITIVE  DIRECTIONS  OF  EXCITATIONS  AND  DISPLACEMENTS 
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FIGURE  2.  (a)  LAMINATED  COMPOSITE  COMPRISING  A  STACK  OF  THIN  FIBRE 
REINFORCED  PLATES;  (b)  STACKING  ANGLE  RELATIVE  TO  LOCAL  (jc,>)  AND 

GLOBAL  (xj)  COORDINATES,  WHICH  ARE  UNRELATED  TO  THOSE  IN  FIGURE  I. 
ACTUAL  CYLINDER  ((^,z)  COORDINATES  ARE  IN  PARENTHESES 
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FIGURE  3.  (a)  CONICAL  SHELL  RNITE  ELEMENT  (WITH  NODES  AT  i  AND  j) 
AND  ITS  DISPLACEMENTS  IN  LOCAL  COORDINATES;  (b)  EXAMPLES  OF  RIB 
CROSS-SECTION  ASSEMBLED  FROM  CONICAL  SHELL  ELEMENTS 
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FIGURE  4.  ACOUSTIC  POWER  IN  WATER  DUE  TO - RADIAL  AND . 

AXIAL  POINT  FORCE  EXCITATION  OF  AIR  FILLED  STEEL  CYLINDER  WITH 
STEEL  T-RIBS;  (a)  NO  RIBS;  (b)  RIBS,  FORCE  ON  SHELL  AT  RIB 
CONNECTION  POINT;  (c)  RIBS,  FORCE  ON  RIB  AT  FLANGE  CENTRE 
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FIGURE  5.  ACOUSTIC  POWER  IN  WATER  DUE  TO - RADIAL  AND . 

AXIAL  POINT  FORCE  EXCITATION  OF  AIR  FILLED  AXIALLY  REINFORCED 
CFRP  CYLINDER  WITH  CIRCUMFERENTIALLY  REINFORCED  CFRP  T-RIBS; 
(a)  NO  RIBS;  (b)  RIBS,  FORCE  AT  RIB  CONNECTION  POINT; 

(c)  RIBS,  FORCE  ON  RIB  AT  FLANGE  CENTRE 
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Annex  A 


SPECTRAL  DYNAMIC  RECEPTANCE  OF  SHELL 
Al.  Differential  Equations 

The  shell  equations  of  motion  are  expressed  in  operator  foim  as 


L,.  Ly  L,j 

E,CM 

0 

Aa  A23 

m 

-  0 

^1  Au  ^ 

PSa,^id-pffl,4>A 

where  E^,  and  are  the  excitation  stresses  applied  to  the  midsurface  of  the  shell;  they  are 
positive  when  acting  in  the  positive  direction  of  the  coordinate  axes.  The  quantities  and  p,  are 
the  pressures  due  to  external  and  internal  fluids,  respectively.  The  operators  Ly  depend  on  the 
particular  shell  theory  used.  Here  the  laminated  composite  shallow  shell  equations  of  Leissa  and  Qatu 
[14],  in  which  effects  due  to  transverse  shear  and  rotary  inertia  are  neglected,  are  adopted.  The 
operators,  in  cylindrical  coordinates  rather  than  in  the  general  coordinates  used  by  these  authors,  are 

L,j  — [  +  TA^^Jadzd^  +  AjjS*/**]  +  p,h^ldt\ 

♦  (A^^^A^)^lad4>dz  *  A^cFldz^], 

I„  +  3B^SPladz‘^d4f  *  *  B^^^ldz*] 

-  ^jja/aSr  ♦  Aj^a’a^j, 

I5,  .-[4„a*/a*a^  +  2A^^^|a^(P^z  >  A^idz^]  * 

Lj^  -+[B„a*/a*a^  ♦  SBj^a’/a’a^az  +  ar*  +  Bj^a’/ar’] 

-  +  A^^a^z^ 

-+[D,ja«/a^a^  +  AD^^Ia^did(f?  +  *  4Djja*/aa^*  + 

*  AJa^-  [2B„a*/a*a^  +  AB^f^|a^^z^^  *  2B^^i^^a^z^  *  P,hcFldt\ 

L21  •  Ljl  •  -L^y  Ijj  ■  -Ljy 

(A.1.2) 
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where  p,  is  the  mean  density  of  the  shell  and  A  is  its  thickness.  When  the  shell  consists  of  a  single 
isotropic  layer  it  can  be  shown  that  these  operators  reduce  to  those  of  the  baseline  cylindrical  shell 
theory  of  Donnell  and  Mushtarl  [8]. 

A2.  Stiffness  Values 

The  stiffness  values  By  and  Dy  are  the  same  as  those  derived  from  laminated  composite  plate 
theory.  The  laminated  plate  comprises  M  thin  layers,  see  Figure  2,  of  fibre  reinforced  orthotropic 
material.  For  the  equivalent  lamination  on  the  shell’s  surface  identify  the  global  x  and  y  axis  with 
the  cylindrical  ^  and  z  axes,  respectively.  In  the  local  x-y  plane  a  typical  layer  has  "engineering* 
elastic  constants:  £„  and  the  Young’s  moduli;  the  shear  modulus;  v^  and  v^  the  Poisson’s 
ratios  which  are  related  by  the  symmetry  relation  For  unidirection  fibre 

reinforcement  in  an  isotropic  matrix  material  simple  mixture  theories  such  as  the  one  proposed  by 
Chamis  [15]  enable  these  constants  to  be  found  from  the  elastic  constants  of  the  component  materials. 

In  the  global  x-y  axes  set  the  orthotropic  constants  of  a  layer  are  obtained  from  the  plane  stress 
matrix  relation 


Oil  Qxx  <?w 

X 

Qn  Qn  Qjx 

t 

Qit  Qtx 

(A.2.1) 


If  the  layer  is  stacked  at  an  angle  $,  the  elements  of  this  constitutive  matrix  are  obtained  from 
standard  transformations,  see,  for  example,  Jones  [16], 


where 


Q^^  -  Q„cob^5  +  2(Q,j+2Q^)sm*5co6*C  + 

Qa  “  +  Ou(sin*C+co««0, 

“  «?u-<?u-2<?«6)sin5cos’C  +  «?u-Oj2+2Q<<)sin*Ccos5, 

(A.2.2) 

Qj,  -  <?„8in"«  *  2«?u+2<?«)fim»Ccoi»C  ♦ 

■  (Cn~Cu“20jj)sin*5cosC  +  (Qi3~Q)3'*'2Q^)sin(co6^(, 


<?„-£,/(l-v^v^), 


are  the  only  non-zero  elements  of  the  constitutive  matrix  in  local  coordinates. 
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The  stiffness  and  density  constants  required  in  equation  (A.  1.2)  are  obtained  by  summing  the  layer 
constants  as  follows 


3a*l 


(A.2.3) 


where  with  z,  and  z^.j  being  the  z  coordinates  of  the  upper  and  lower  surfaces  of  the 

(m)th  layer;  z,  is  the  z  coordinate  of  tlie  laminate's  midsurface;  p  is  the  layer  density;  h  is  the  total 
thickness  of  the  laminate.  The  stiffness  values  are  related  to  in-plane  stretching  of  the  shell's 
midsurface;  the  values  to  bending  of  the  surface;  the  values  to  coupling  between  bending  and 
stretching.  When  the  shell  is  laminated  symmetrically  with  respect  to  its  midsurface,  the  values 
are  zero. 


A3.  Spectral  Dynamic  Stiffness 

Rqplace  the  field  quantities  in  equation  (A.  1.1)  by  their  Fourier  transforms,  defined  by  equation 
(2.1.1),  to  give  the  spectral  dynamic  stiffness  matrix  relation 


Sj,(n,«)  Sjj(n,«)  S^(H,a) 
Sjj(n,a)  Sjj(fi,a)  S„(n,a) 


«/»,«) 


BJifha) 

£^(n,o) 


(A.3.1) 


where 


♦  [AgjiiVa*  lA^ernfa  *  Ajja*]  -  p/«*, 

S,,  -  ♦  ♦  (Au+AJo/i/fl  ♦ 

Sj,  ■  -f  ♦  3£jjaVa  +  +  £jjo’  ♦  Ay^ala  +  i4,,»t/a*],(A-3.2) 

Sj2  ■  +  *  lAy^tila  *  -  p/«’, 

Sjj  ■  -f  ♦  3£,j»*o/a’  +  (B^.  *  IBf^a^nfa  r  Bj^a}  ♦  A„n/o* 
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5jj  ■  [DjjnVa^  +  4D,jn*o/a*  ♦  2(Djj  +  +  AD^^na^la 

*  D^a"  +  2B^^n^|a^  +  45,jna/fl*  +  IByjO^la  +  -  p/iu* 

^21  ■  ^U*  ^31  *  “^13*  ^32  “  "^23* 

The  last  two  terms  of  are  the  fluid  loading  spectral  coefficients  of  the  exterior  and  interior  fluids, 
respectively,  see,  for  example,  reference  [7]. 

The  solution  of  equation  (A.3.1)  is  the  spectral  receptance  matrix  relating  spectral  forces  and 
displacements,  viz. 


Djj(rt,a)  Djj(«,a)  D^^(n,ol) 

EJinfi) 

«/».«) 

Z>j,(n,a)  Djj(n,a)  D^(n,a) 

2)jj(n,o)  D^^(n,a)  D„(n,a) 

EJinfl) 

As  a  rib  is  assumed  to  exert  a  meridional  moment  on  the  shell,  in  addition  to  the  three  orthogonal 
stresses,  it  is  necessary  to  expand  equation  (A.3.3)  to  include  moment  excitation  per  unit  area, 
via  its  spectral  quantity,  Af^(n,a) .  The  moment  excitation  is  equivalent  to  a  radial  stress 
excitation  -9M^(0,z)/az,  whose  spectral  form  is  -iaAf^(n,a);  the  corresponding  meridional  rotation, 
has  spectral  form  T(n,a)»io«,(n,o).  Thus,  noting  that  the  elements  2)^  are 

interpreted  as  the  spectral  response  in  direction  i  due  to  unit  spectral  excitation  in  direction/,  it  can 
be  shown  that  the  expanded  matrix  relation  is 


«/n,o) 

D^^(n,a)  2Jjj(n,a)  D,j(n,o)  D^^in,a) 

E^(n,a) 

2)j,(rt,a)  l>2j(n,o)  Djj(n,a) 

E/n.a) 

■ 

DJ^{n,a)  2?3,(n,a)  Dy^{n,a.) 

E/n.a) 

D^^(n.a)  D^(n,o)  D^j(n,a)  2>^(n,a) 

A//n,a) 

A^ 


where 


D^/,n,a)  -  -iaD^^(n,a), 
D^(na)  -  -ioZ)jj(n,a), 
I>4a(».«)  -  ittX>ji(n,«). 
i>44(«»<‘)  “  a*i)„(/«,a). 


024(«.«)  "  -iaD^itua), 
^4i(n.«)  “  ia2)j,(n.a). 


Denote  the  4  x  4  coefficient  matrix  ty  D(/i.o) :  this  is  the  shell’s  spectral  dynamic  receptance  matrix 
which  is  required  in  equation  (2.2.1). 
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Annex  B 


DYNAMIC  STIFFNESS  MATRK  OF  RIB 
Bl.  Finite  Element  Modelling  of  Rib 

A  rib  is  modelled  as  an  assemblage  of  axisymmetric  conical  shell  elements,  whose  versatility  is 
illustrated  in  Figure  3  which  shows  examples  of  rib  cross-section  which  can  be  generated  from  these 
elements. 

The  theory  of  a  conical  shell  element  is  based  on  the  strain  definitions  in  the  local  surface  coordinate 
system  (^,s) 

•  du^frd^  *  iujxa^  ♦  H^sinP)/r, 
e,  -  dujds, 

*  du^ds  -  u^inp/r, 

(B.1.1) 

+  cospa<^/r^5^  -  ^^dujrds, 

Xt  • 

t  «  -l^ujrdsd^  +  2sinP(3iu,/r*30  *  2cos.^(du^lrds)  -  2sinPcosP(u^/r*), 

I 

where  and  are  raidsurface  normal  strains;  is  the  shear  strain;  and  x,  are  cha:.ges  in 
curvature;  t  is  twist;  and  are  displacements  normal  to  the  surface,  around  the 

circumference  and  along  the  generator,  respectively;  p  is  the  slope  of  the  conical  shell  generator. 

j 

The  6x6  constitutive  matrix  D^,  relating  stress  and  moment  resultants  to  midsurface  strain,  is 
defined  by  the  matrix  relation 


•^11  '^l*  ^11 

A],  Ajf  ^  Bjs 

^  ^«6 

y.* 

^11  ^12 

K 

^  ^3*  ^32  ^3* 

X, 

^3t 

T 

where  the  stiffness  coefficients  have  been  defmed  in  Annex  A. 
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B2.  Mass  and  Stiffness  Matrices  of  Element 

Field  quantities  are  expanded  as  Fourier  series  in  the  circumferential  coordinate.  Thus,  for  the 
function 


m 

T 


where  POi^)  -  (l/2n)  . 


(B.2.1) 


Because  of  the  orthogonality  property  of  the  harmonic  functions  a  finite  element  analysis  can  be 
carried  out  for  each  harmonic  P(n^)  separately.  Thus,  the  following  procedure  refers  to  values  in 
the  nth  harmonic,  and  the  requirement  to  sum  the  Fourier  series,  over  values  of  n,  to  obtain  tiie  actual 
field  quantity  is  implicitly  assumed.  The  theory  is  given  with  minimal  explanation:  details  for  cos(n<^) 
variation  are  contained  in  the  pioneering  p^er  of  Percy  et  al.  [17]  and  numerous  texts  on  the  finite 
element  method;  extension  of  exp(in^)  variation  as  herein  is  not  difficult. 


The  conical  shell  displacement  uj[n^)  normal  to  its  generator  is  represented  by  a  cubic  in  the  variable 
s  which  is  the  distance  along  the  generator  from  node  i .  The  displacement  along  the  generator 
ujin^)  and  the  displacement  around  the  circumference  u^(n,s)  are  both  assumed  to  vary  linearly  in 
the  coordinate  s.  The  rotation  of  the  meridian  is  W .  The  equation  of  motion  of 
the  element  is  found  by  standard  finite  element  procedure  as 


where 


S,(n)-o)*M,]u,(n)  -  E,(n), 


S,(n)  -  2n(.\;y  [  /^B>A)D,B,(n,j)r&]A;‘, 
M.  .  27t(A;V  [/^^pAV,Wr&]A;'. 


(B.2.2) 


(B.2.3) 


The  8x8  matrix  S,(/i)  is  the  stiffness  matrix  of  the  element  in  the  nth  harmonic;  the  8  x  8  matrix  M, 
is  the  mass  matrix  which  is  independent  of  the  harmonic  number;  the  superscripts  -1,  *  and  t  denote 
the  operations  of  inverse,  complex  conjugate  and  transpose,  respectively;  the  subscript  e  indicates  that 
the  quantity  to  which  it  is  attached  refers  to  a  single  conical  shell  element;  the  matrices  B, ,  and 
are  given  below;  the  6  x  6  matrix  D,  is  the  constitutive  matrix  of  equation  (B.1.2);  p,  is  the 
element’s  density  and  \  is  its  thickness.  The  integrals  are  evaluated  numerically,  by  a  simple 
quadrature  scheme,  to  avoid  lengthy  algebra. 


The  matrices  u/n)  and  E,(n)  are  8  x  1  column  vectors  of  nodal  displacements  and  nodal  excitations 
in  the  global  coordinate  system,  viz 


E,(n)  -  {<(n),£>),£;(«),*f>),£/(n),£;(n),£/(n^ 
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where  the  superscripts  i  and  j  refer  to  the  left  and  right  nodes  on  the  element.  For  a  force  vector, 
>  per  unit  length  in  the  circumferential  direction  at  node  i ,  the  excitation  matrix  at  this  node 
can  be  shown  to  be  E'(n)-2itr^^n),  where  F/n)  is  obtained  from  the  Fourier  transform  equation 
(B.2.1)  and  is  the  radial  coordinate  of  the  node;  the  excitation  matrix  E^(n)  at  node  /  is  identically 
zero  if  there  is  no  excitation  there.  Thus,  for  a  point  force  applied  at  node  i  at 

the  4  X  1  excitation  matrix  is  simply  E'(n)-Fexp(-in^^  where  the  matrix  F  is  defined  in 
Section  2.4. 

The  non-zero  elements  of  the  6  x  8  matrix  B/n)  are 

(l/r)sinp,  5^=  (j/r)smp,  in/r,  Uts[r, 

cosp/r,  B,j=  (s/r)cosp,  (5Vr)coep,  (?Vr)cosp, 

Bjj-  1,  Bjj*  Utir,  Bjj-  insir,  -sinP/r,  l-(j/r)smp, 

(in/f*)co«P,  B^^~  (ins/r*)cosP,  nVr*. 

B44-  (fi*s/r‘*)-sinP/r,  (n^sV*)~(24)sinP/r,  (B.2.5) 

^41*  -2  Bj,»  -d«, 

-(2/r*)sinPcosp,  (2/r)ccsP-(2$/r*)sinPcosP, 

Bj,»  2iRsiaP/r^  Bj4=  -(2in/r)+(2ins/r*)sinP, 

B„"  -(4i/«s/r)+(2insVr*)sinP,  Bj,=  -(6/nsVr)+(2insV'’*)sinP, 
and  the  non-zero  elements  of  the  8  x  8  matrix  are 

A||-  cosp,  -sinp,  1,  A,,-  sinp,  A„>  cosp, 

A^j-  1,  Ajj*  cosp,  Ajj-  IcosP,  Aj5*  -amp,  Aj^*  -IsinP,  (B.2.6) 

A„«  -L^P,  Aj|«  -L’sinP,  A^,*  1,  A^»  L,  A,,»  sinP, 
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IsinP,  cosP,  Ajf=  Lcosfi,  L^coefi,  A^a  L’cosP, 
A„a  1.  A„=  21,  A„=  3L\ 


while  the  non-zero  elements  of  the  8  x  8  matrix  V^,  rotary  inertia  being  ignored,  are 


*'u=S. 

K„=l. 

1,  s,  s*,  Kj,®  5®, 

5*,  K^=  s*,  K,j=  s\  s\ 

y„.s’,  ('.-s*. 


(B.2.7) 


B3.  Dynamic  Stiffness  of  Rib 

Assembly  of  a  number  of  conical  shell  elements  to  form  a  rib,  as  shown  in  Figure  3,  for  example, 
is  a  straightforward  finite  element  computational  procedure  which  reflects  continuity  of  displacement 
and  equilibrium  of  stresses  at  the  nodes.  When  the  rib  has  N finite  element  nodes  the  assembly  results 
in  the  system  matrix  equation 


Z,(n)U,(«)  -  E,(n), 


(B.3.1) 


where  Z/n)  is  a  4Wx  4^  system  dynamic  stiffness  matrix,  U,.(n)  is  a  4Af  x  1  column  containing  the 
displacements  of  all  nodes  and  E/n)  is  a  4iV  x  1  column  vector  containing  the  nodal  excitations.  The 
subscript  r  refers  to  a  complete  rib.  Take  the  inverse  of  equation  (B.3.1)  to  give  the  system  dynamic 
receptance  matrix  Z/n)*'  whose  ijtb  element  is  interpreted  as  the  response  at  degree  of  freedom  i  due 
to  unit  excitation  at  degree  of  freedom  J.  The  inverse  can  usually  be  computed  rapidly  if  the  elements 
of  Z^(n)  are  stored  in  a  format  suitable  for  a  band  solver. 

If  a  rib  is  attached  to  the  shell  at  node  p  (usually  1)  and  external  excitation  forces  are  to  be  applied 
at  node  g,  then  delete  all  rows  and  columns  of  Z,(n)~‘  not  associated  with  the  eight  degrees  of 
freedom  at  nodes  p  and  9.  The  following  partitioned  recqitance  matrix  results 


V") 

V") 

E/n) 

n,(») 

R,(«)  R„(") 

E,(») 

where  u^(fi)  and  u^(n)  are  4  x  1  column  vectors  of  dbplacements  at  nodes  p  and  respectively; 
E^(fi)  and  E^(n)  are  column  vectors  of  internal  reaction  forces  and  external  excitations;  R^(n)  is  a 
4x4  recqitance  matrix,  whose  columns  are  the  4  x  1  response  vector  at  p  due  to  unit  excitation  in 
turn  at  each  of  the  four  degrees  of  freedom  at  q. 


B-4 


Equation  (B.3.2)  can  be  written  as  the  two  matrix  relations 


u/«)  -  R^(«)E/n)  ♦  R^(ii)E,(n). 

u,(n)  -  R,(n)E/n)  *  R^(n)E,(n) 

The  first  of  these  equations  is  re-arranged  as 

F/ii)  -  (l/2xrpR^(n)-'u//i)  -  (^iir/Lnr^V^(nr*R^(n)F^(n), 


(B.3.3) 


(B.3.4) 


for  line  force  excitations  P^(^)  and  F^(4>)  at  nodes  p  and  q,  respectively,  F^(n)  and  F^(n)  being 
the  Fourier  series  transforms  according  to  equation  (B.2.1). 

In  the  absence  of  external  excitation  on  a  rib,  F^(n)-0,  the  dynamic  stiffness  matrix  of  the  rib  has 
been  defined  in  [3]  as  F^(n)-B(n)u^(n).  Noting  that  the  rib  attachment  point  is  r^*a,  this  gives 
B(n)»(l/2iia)R^(«)‘*,  which  matrix  is  required  in  the  shell’s  spectral  response  of  equations  (2.1.1), 

The  second  term  in  (B.3.4)  is  that  part  of  the  reaction  force  at  point  p  which  is  due  to  line  force 
excitation  F^(n)  on  the  same  rib.  The  equivalent  excitation  on  the  shell’s  surface  must  be  equal  but 
of  opposite  sign.  Thus 

E(n,z)  -  (2«r;2nn)R^(/i)-'R^(n)F,(n)8(z-z^.  (B.3.5) 

where  the  rib  attachment  point  is  at  r-o,  z-Zg.  The  delta  function  has  dimension  m*'  in  order  to 
convert  forces  per  unit  length  in  the  circumferential  direction  on  the  rib  into  the  required  forces  per 
unit  area  on  the  shell.  Take  the  Fourier  integral  transform  of  (B.3.S)  and  let  the  excitation  be  the 
point  force  vector  F^(^)-F5(^-^^/r^  whose  Fourier  series  transform  is  F^(n)-Fexp(-in^^/2icr^, 
to  give 


E(n,a)  -  R^(n)*'R^(/i)Fexp(-in0,-ioz^/2«a,  (B.3.6) 

which  is  the  equivalent  shell  excitation  required  in  Section  2.S. 

The  second  of  equations  (B.3.3)  can  be  used  to  calculate  the  displacement  Uy(n)  at  the  excitation  point 
on  a  rib,  but  this  is  a  formidable  problem  as  it  requires  the  displacement  of  the  shell  u^(n,z-z^  at  the 
rib  attachmeU  point,  which  can  only  be  obtained  by  evaluating  a  Fourier  integral  numerically. 


i 
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